Computational Structure Prediction Methods Research Articles

Current Computational Methods for Protein-Peptide Complex Structure Prediction.

Peptide-mediated protein-protein interactions (PPIs) play an important role in various biological processes. The development of peptide-based drugs to modulate PPIs has attracted increasing attention due to the advantages of high specificity and low toxicity. In the development of peptide-based drugs, one of the most important steps is to determine the interaction details between the peptide and the target protein. In addition to experimental methods, recently developed computational methods provide a cost-effective way for studying protein-peptide interactions. In this article, we carefully reviewed recently developed protein-peptide docking methods, which were classified into three groups: template-based docking, template-free docking, and hybrid method. Then, we presented available benchmarking sets and evaluation metrics for assessing protein-peptide docking performance. Furthermore, we discussed the use of molecular dynamics simulations, as well as deep learning approaches in protein-peptide complex prediction.

Protein Secondary Structure Prediction Using Cascaded Feature Learning Model

The protein secondary structure prediction (PSSP) is pivotal for predicting tertiary structure, which is proliferating in demand for drug design and development. Further, it can be used to learn different protein functions. Although there are many computational methods for protein structure prediction, none of them have succeeded 100% in solving the protein structure prediction problem. This paper proposes a deep learning model, namely Cascaded Feature Learning Model (CFLM) that uses a multi-stage transfer learning approach based on the Residual Dense Network (RDN) for predicting protein secondary structure. The model is trained with different protein datasets at various levels of transfer learning and accepts selected protein features such as solvent accessibility, physicochemical properties, PSSM (position-specific scoring matrix), and PSFM (position-specific frequency matrix) as input. The Q3 and Q8 accuracy obtained on CASP 13 and 14 benchmark datasets with the proposed approach are 91.23% and 91.45%, and 76.83% and 78.04%, respectively. The performance of CFLM is also compared with some of the recent PSSP approaches and it is observed that the proposed approach improves the prediction accuracy by 2.8% in terms of Q3 and 1.81% in terms of Q8 in the CASP 14 dataset.

Pseudoknots in RNA Structure Prediction.

RNA molecules play active roles in the cell and are important for numerous applications in biotechnology and medicine. The function of an RNA molecule stems from its structure. RNA structure determination is time consuming, challenging, and expensive using experimental methods. Thus, much research has been directed at RNA structure prediction through computational means. Many of these methods focus primarily on the secondary structure of the molecule, ignoring the possibility of pseudoknotted structures. However, pseudoknots are known to play functional roles in many RNA molecules or in their method of interaction with other molecules. Improving the accuracy and efficiency of computational methods that predict pseudoknots is an ongoing challenge for single RNA molecules, RNA-RNA interactions, and RNA-protein interactions. To improve the accuracy of prediction, many methods focus on specific applications while restricting the length and the class of the pseudoknotted structures they can identify. In recent years, computational methods for structure prediction have begun to catch up with the impressive developments seen in biotechnology. Here, we provide a non-comprehensive overview of available pseudoknot prediction methods and their best-use cases. © 2023 Wiley Periodicals LLC.

Fast Local Alignment of Protein Pockets (FLAPP): A System-Compiled Program for Large-Scale Binding Site Alignment

Protein function is a direct consequence of its sequence, structure, and the arrangement at the binding site. Bioinformatics using sequence analysis is typically used to gain a first insight into protein function. Protein structures, on the other hand, provide a higher resolution platform into understanding functions. As the protein structural information is increasingly becoming available through experimental structure determination and through advances in computational methods for structure prediction, the opportunity to utilize these data is also increasing. Structural analysis of small molecule ligand binding sites in particular provides a direct and more accurate window to infer protein function. However, it remains a poorly utilized resource due to the huge computational cost of existing methods that make large-scale structural comparisons of binding sites prohibitive. Here, we present an algorithm called FLAPP that produces very rapid atomic level alignments. By combining clique matching in graphs and the power of modern CPU architectures, FLAPP aligns a typical pair of binding sites at ∼12.5 ms using a single CPU core, ∼1 ms using 12 cores on a standard desktop machine, and performs a PDB-wide scan in 1-2 min. We perform rigorous validation of the algorithm at multiple levels of complexity and show that FLAPP provides accurate alignments. We also present a case study involving vitamin B12 binding sites to showcase the usefulness of FLAPP for performing an exhaustive alignment-based PDB-wide scan. We expect that this tool will be invaluable to the scientific community to quickly align millions of site pairs on a normal desktop machine to gain insights into protein function and drug discovery for drug target and off-target identification and polypharmacology.

Unravelling the Origin of Solvate Formation in the Anticancer Drug Trametinib: Insights from Crystal Structure Analysis and Computational Modeling

The discovery of new crystalline forms of large, flexible drug molecules is of significant interest to the pharmaceutical industry as a means of improving one or more physicochemical properties of the active pharmaceutical ingredient (API). Trametinib (TRMB) is an anticancer drug, for which no solid-state forms are yet reported. Here, we report the results of an extensive screen for cocrystals of TRMB in an effort to improve the poor aqueous solubility of this API. Contrary to expectations, this effort led to the discovery of only six isostructural solvates of TRMB, where the solvent molecules fill the voids of the crystal structure and are stabilized by various intermolecular interactions. Computational crystal structure prediction methods were used to rationalize the high propensity for solvent inclusion in TRMB as arising from the poor crystal packing in TRMB. All predicted polymorphs of TRMB within a relative lattice energy of 44 kJ mol–1 of the global minimum display solvent accessible voids ranging from 1.4 to 16.5% of the unit cell volume. The serendipitous discovery of solvates in experiments targeting cocrystals,strongly suggests that the well-established synthon approach for the discovery of cocrystals works relatively well for documented literature examples that use rigid coformers with minimal conformational flexibility. However, for large APIs such as TRMB, in-silico screening methods should be adopted to reduce the need for extensive experimental cocrystal screens as the ability to form a cocrystal for drug candidates relies on more than synthon complementarity.

Computational Pipeline for Reference-Free Comparative Analysis of RNA 3D Structures Applied to SARS-CoV-2 UTR Models.

RNA is a unique biomolecule that is involved in a variety of fundamental biological functions, all of which depend solely on its structure and dynamics. Since the experimental determination of crystal RNA structures is laborious, computational 3D structure prediction methods are experiencing an ongoing and thriving development. Such methods can lead to many models; thus, it is necessary to build comparisons and extract common structural motifs for further medical or biological studies. Here, we introduce a computational pipeline dedicated to reference-free high-throughput comparative analysis of 3D RNA structures. We show its application in the RNA-Puzzles challenge, in which five participating groups attempted to predict the three-dimensional structures of 5- and 3-untranslated regions (UTRs) of the SARS-CoV-2 genome. We report the results of this puzzle and discuss the structural motifs obtained from the analysis. All simulated models and tools incorporated into the pipeline are open to scientific and academic use.

The integration of AlphaFold-predicted and crystal structures of human trans-3-hydroxy-l-proline dehydratase reveals a regulatory catalytic mechanism

Computational methods for protein structure prediction have made significant strides forward, as evidenced by the last development of the neural network AlphaFold, which outperformed the CASP14 competitors by consistently predicting the structure of target proteins. Here we show an integrated structural investigation that combines the AlphaFold and crystal structures of human trans-3-Hydroxy-l-proline dehydratase, an enzyme involved in hydroxyproline catabolism and whose structure had never been reported before, identifying a structural element, absent in the AlphaFold model but present in the crystal structure, that was subsequently proved to be functionally relevant. Although the AlphaFold model lacked information on protein oligomerization, the native dimer was reconstructed using template-based and ab initio computational approaches. Moreover, molecular phasing of the diffraction data using the AlphaFold model resulted in dimer reconstruction and straightforward structure solution. Our work adds to the integration of AlphaFold with experimental structural and functional data for protein analysis, crystallographic phasing and structure solution.

Solid Nitrogen and Nitrogen‐Rich Compounds as High‐Energy‐Density Materials

Polymeric allotropes of nitrogen and nitrogen‐rich compounds containing NN single bonds are attractive candidates for high‐energy‐density materials. However, the low kinetic stability of single bond in nitrogen means that these materials have to be synthesized under high pressure and temperature. Recent advances of experimental high‐pressure techniques and computational structure prediction methods have enabled the investigation of many new nitrogen phases and nitrogen‐rich compounds, and some have been successfully synthesized. Some of the recently synthesized materials have been retrieved at ambient conditions, and the large hysteresis warrants further application‐inspired research. This review introduces the important classes of nitrogen‐rich materials ranging from fundamental insight into bonding in nitrogen, to theoretical investigation of hypothetical compounds, and to laboratory realization of these energetic materials in diamond anvil cell. The progressive development of the field is narrated in chronological order, with the emphasis placed on recent progress in the last 2 years.

Tuning intrinsic disorder predictors for virus proteins.

Many virus-encoded proteins have intrinsically disordered regions that lack a stable, folded three-dimensional structure. These disordered proteins often play important functional roles in virus replication, such as down-regulating host defense mechanisms. With the widespread availability of next-generation sequencing, the number of new virus genomes with predicted open reading frames is rapidly outpacing our capacity for directly characterizing protein structures through crystallography. Hence, computational methods for structural prediction play an important role. A large number of predictors focus on the problem of classifying residues into ordered and disordered regions, and these methods tend to be validated on a diverse training set of proteins from eukaryotes, prokaryotes, and viruses. In this study, we investigate whether some predictors outperform others in the context of virus proteins and compared our findings with data from non-viral proteins. We evaluate the prediction accuracy of 21 methods, many of which are only available as web applications, on a curated set of 126 proteins encoded by viruses. Furthermore, we apply a random forest classifier to these predictor outputs. Based on cross-validation experiments, this ensemble approach confers a substantial improvement in accuracy, e.g., a mean 36 per cent gain in Matthews correlation coefficient. Lastly, we apply the random forest predictor to severe acute respiratory syndrome coronavirus 2 ORF6, an accessory gene that encodes a short (61 AA) and moderately disordered protein that inhibits the host innate immune response. We show that disorder prediction methods perform differently for viral and non-viral proteins, and that an ensemble approach can yield more robust and accurate predictions.

Carbohydrate-Protein Interactions: Advances and Challenges.

A carbohydrate, also called saccharide in biochemistry, is a biomolecule consisting of carbon (C), hydrogen (H) and oxygen (O) atoms. For example, sugars are low molecular-weight carbohydrates, and starches are high molecular-weight carbohydrates. Carbohydrates are the most abundant organic substances in nature and essential constituents of all living things. Protein-carbohydrate interactions play important roles in many biological processes, such as cell growth, differentiation, and aggregation. They also have broad applications in pharmaceutical drug design. In this review, we will summarize the characteristic features of protein-carbohydrate interactions and review the computational methods for structure prediction, energy calculations, and kinetic studies of protein-carbohydrate complexes. Finally, we will discuss the challenges in this field.

Hierarchical, rotation-equivariant neural networks to select structural models of protein complexes.

Predicting the structure of multi-protein complexes is a grand challenge in biochemistry, with major implications for basic science and drug discovery. Computational structure prediction methods generally leverage predefined structural features to distinguish accurate structural models from less accurate ones. This raises the question of whether it is possible to learn characteristics of accurate models directly from atomic coordinates of protein complexes, with no prior assumptions. Here we introduce a machine learning method that learns directly from the 3D positions of all atoms to identify accurate models of protein complexes, without using any precomputed physics-inspired or statistical terms. Our neural network architecture combines multiple ingredients that together enable end-to-end learning from molecular structures containing tens of thousands of atoms: a point-based representation of atoms, equivariance with respect to rotation and translation, local convolutions, and hierarchical subsampling operations. When used in combination with previously developed scoring functions, our network substantially improves the identification of accurate structural models among a large set of possible models. Our network can also be used to predict the accuracy of a given structural model in absolute terms. The architecture we present is readily applicable to other tasks involving learning on 3D structures of large atomic systems.

Hybrid methods for combined experimental and computational determination of protein structure.

Knowledge of protein structure is paramount to the understanding of biological function, developing new therapeutics, and making detailed mechanistic hypotheses. Therefore, methods to accurately elucidate three-dimensional structures of proteins are in high demand. While there are a few experimental techniques that can routinely provide high-resolution structures, such as x-ray crystallography, nuclear magnetic resonance (NMR), and cryo-EM, which have been developed to determine the structures of proteins, these techniques each have shortcomings and thus cannot be used in all cases. However, additionally, a large number of experimental techniques that provide some structural information, but not enough to assign atomic positions with high certainty have been developed. These methods offer sparse experimental data, which can also be noisy and inaccurate in some instances. In cases where it is not possible to determine the structure of a protein experimentally, computational structure prediction methods can be used as an alternative. Although computational methods can be performed without any experimental data in a large number of studies, inclusion of sparse experimental data into these prediction methods has yielded significant improvement. In this Perspective, we cover many of the successes of integrative modeling, computational modeling with experimental data, specifically for protein folding, protein-protein docking, and molecular dynamics simulations. We describe methods that incorporate sparse data from cryo-EM, NMR, mass spectrometry, electron paramagnetic resonance, small-angle x-ray scattering, Förster resonance energy transfer, and genetic sequence covariation. Finally, we highlight some of the major challenges in the field as well as possible future directions.

Computational insight into the three-dimensional structure of ADP ribosylation factor like protein 15, a novel susceptibility gene for rheumatoid arthritis

The ARL15 gene (ADP ribosylation factor like protein 15) encodes for an uncharacterized small GTP-binding protein. Its exact role in human physiology remains unknown, but a number of genetic association studies have recognised different variants in this gene to be statistically associated with numerous traits and complex diseases. We have previously reported a novel association of ARL15 with rheumatoid arthritis (RA) based on a genome-wide association study in a north Indian cohort. Subsequent investigations have provided leads for its involvement in RA pathophysiology, especially its potential as a novel therapeutic target. However, the absence of an experimentally determined tertiary structure for ARL15 significantly hinders the understanding of its biochemical and physiological functions, as well as development of potential lead molecules. We, therefore, aimed to derive a high quality, refined model of the three dimensional structure of human ARL15 protein using two different computational protein structure prediction methods – template-based threading and ab initio modelling. The best model each from among the five each derived from both the approaches was selected based on stringent quality assessment and refinement. Molecular dynamics simulations over long timescales revealed the ab initio model to be relatively more stable, and it marginally outperformed the template-based model in the quality assessment as well. A putative GTP-binding site was also predicted using homology for the ARL15 protein, where potential competitive inhibitors can be targeted. This high quality predicted model may provide insights to the biological role(s) of ARL15 and inform and guide further experimental, structural and biochemical characterization efforts. Communicated by Ramaswamy H. Sarma

Ab Initio Prediction of Metal-Organic Framework Structures

Metal-organic frameworks (MOFs) have emerged as highly versatile materials with applications in gas storage and separation, solar light energy harvesting and photocatalysis. The design of new MOFs, however, has been hampered by the lack of computational methods for ab initio MOF structure prediction, which could be used to inspire and direct experimental synthesis. Here, we report the first ab initio method for the prediction of MOF structures and test it against a diverse set of known MOFs that were chosen for their differences in topology, metal coordination geometry, and ligand binding sites. In all cases, our calculations produced structures that match experiment using only the target composition and ligand molecular structure, proving the versatility of our procedure. The herein presented methodology utilizes the point group symmetry of ligands to enable, for the first time, prediction of MOF structures from first principles, without having to resort to empirical guidelines based on rigid connectivity of nodes and linkers, or to previously determined crystal structures and topologies of known MOFs. This advance provides the first tool to change MOF design from an empirically based process that is based on chemistʼs intuition rooted in literature- or database-established knowledge of node-and-linker connectivity to a more general and theory-driven materials development. This ab initio MOF structure prediction approach, which is here validated on a range of known MOF classes, provides a unique opportunity to explore the phase landscape of MOFs computationally and enables MOF research and development even in case of limited access to laboratory resources, as for example in case of a global pandemic.

Multimedia Jupyter Notebooks for Learning Structure Prediction and Design

PyRosetta is a Python library of the Rosetta computational methods for protein structure prediction and design. Because the field attracts students from many different majors, teaching a class with a wide variety of skills and experience levels is a challenge. Therefore, the goal of this project is to develop interactive multimedia material that promotes self-paced active learning. We have developed a hands-on education strategy with a set of nine modules to teach topics in the field, from protein structural analysis to protein folding and docking. These workshop modules are created in Jupyter Notebooks, a shareable web application that supports live code, visualizations, and text. When used in the Google Colaboratory environment, students can access and learn PyRosetta with no local machine setup necessary. By capitalizing on the advantages of this online digital format, we have embedded images, molecular visualization movies, and interactive coding examples. This multimedia approach may better reach students from different majors and experience levels as well as attract more researchers from smaller labs and related disciplines to leverage PyRosetta in their work.

Predicting Protein Complex Structure from Surface-Induced Dissociation Mass Spectrometry Data.

Recently, mass spectrometry (MS) has become a viable method for elucidation of protein structure. Surface-induced dissociation (SID), colliding multiply charged protein complexes or other ions with a surface, has been paired with native MS to provide useful structural information such as connectivity and topology for many different protein complexes. We recently showed that SID gives information not only on connectivity and topology but also on relative interface strengths. However, SID has not yet been coupled with computational structure prediction methods that could use the sparse information from SID to improve the prediction of quaternary structures, i.e., how protein subunits interact with each other to form complexes. Protein–protein docking, a computational method to predict the quaternary structure of protein complexes, can be used in combination with subunit structures from X-ray crystallography and NMR in situations where it is difficult to obtain an experimental structure of an entire complex. While de novo structure prediction can be successful, many studies have shown that inclusion of experimental data can greatly increase prediction accuracy. In this study, we show that the appearance energy (AE, defined as 10% fragmentation) extracted from SID can be used in combination with Rosetta to successfully evaluate protein–protein docking poses. We developed an improved model to predict measured SID AEs and incorporated this model into a scoring function that combines the RosettaDock scoring function with a novel SID scoring term, which quantifies agreement between experiments and structures generated from RosettaDock. As a proof of principle, we tested the effectiveness of these restraints on 57 systems using ideal SID AE data (AE determined from crystal structures using the predictive model). When theoretical AEs were used, the RMSD of the selected structure improved or stayed the same in 95% of cases. When experimental SID data were incorporated on a different set of systems, the method predicted near-native structures (less than 2 Å root-mean-square deviation, RMSD, from native) for 6/9 tested cases, while unrestrained RosettaDock (without SID data) only predicted 3/9 such cases. Score versus RMSD funnel profiles were also improved when SID data were included. Additionally, we developed a confidence measure to evaluate predicted model quality in the absence of a crystal structure.

Homology modeling in drug discovery: Overview, current applications, and future perspectives.

Homology modeling is one of the computational structure prediction methods that are used to determine protein 3D structure from its amino acid sequence. It is considered to be the most accurate of the computational structure prediction methods. It consists of multiple steps that are straightforward and easy to apply. There are many tools and servers that are used for homology modeling. There is no single modeling program or server which is superior in every aspect to others. Since the functionality of the model depends on the quality of the generated protein 3D structure, maximizing the quality of homology modeling is crucial. Homology modeling has many applications in the drug discovery process. Since drugs interact with receptors that consist mainly of proteins, protein 3D structure determination, and thus homology modeling is important in drug discovery. Accordingly, there has been the clarification of protein interactions using 3D structures of proteins that are built with homology modeling. This contributes to the identification of novel drug candidates. Homology modeling plays an important role in making drug discovery faster, easier, cheaper, and more practical. As new modeling methods and combinations are introduced, the scope of its applications widens.

Exploring the conformational space for protein folding with sequential Monte Carlo

Computational methods for protein structure prediction from amino acid sequence are of vital importance in modern applications, for example protein design in biomedicine. Efficient sampling of conformations according to a given energy function remains a bottleneck, yet is a vital step for energy-based structure prediction methods. While the Protein Data Bank of experimentally determined 3-D protein structures has steadily increased in size, structure predictions for new proteins tend to be unreliable in the amino acid segments where there is low sequence similarity with known structures. In this paper we introduce a new method for building such segments of protein structures, inspired by sequential Monte Carlo methods. We apply our method to examples of real 3-D structure predictions and demonstrate its promise for improving low confidence segments. We also provide applications to the prediction of reconstructed segments in known structures, and to the assessment of energy function accuracy. We find that our method is able to produce conformations that have both low energies and good coverage of the conformational space and hence can be a useful tool for protein design and structure prediction.

Bioinformatics Tools and Benchmarks for Computational Docking and 3D Structure Prediction of RNA-Protein Complexes.

RNA-protein (RNP) interactions play essential roles in many biological processes, such as regulation of co-transcriptional and post-transcriptional gene expression, RNA splicing, transport, storage and stabilization, as well as protein synthesis. An increasing number of RNP structures would aid in a better understanding of these processes. However, due to the technical difficulties associated with experimental determination of macromolecular structures by high-resolution methods, studies on RNP recognition and complex formation present significant challenges. As an alternative, computational prediction of RNP interactions can be carried out. Structural models obtained by theoretical predictive methods are, in general, less reliable compared to models based on experimental measurements but they can be sufficiently accurate to be used as a basis for to formulating functional hypotheses. In this article, we present an overview of computational methods for 3D structure prediction of RNP complexes. We discuss currently available methods for macromolecular docking and for scoring 3D structural models of RNP complexes in particular. Additionally, we also review benchmarks that have been developed to assess the accuracy of these methods.

A method to improve prediction of secondary structure for large single RNA sequences

The function of a particular RNA molecule within an organic system is principally determined by its structure. The current physical methods available for structure determination are time consuming and expensive. Hence, computational methods for structure prediction are often used.The prediction of the structure of a large single sequence of RNA needs a lot of research work. In the present work, a method is introduced to improve the prediction of large single sequence RNA secondary structure obtained by Mfold program using the concept of minimum free energy. The Mfold program contains a constraint option that allows forcing some helices in the predicted structure. In our method, some of the firstly formed hairpins that are expected, by a statistical study, to be present in the real structure are forced in the Mfold predicted structure. The results show improvement, toward the real structure, in the Mfold predicted structure and this gives evidence to the RNA kinetic folding.